Quantum coherence of double-well BEC: a SU(2)-coherent-state path-integral 

approach 
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Macroscopic quantum coherence of Bose gas in a double-well potential is studied based on SU(2)- 
coherent-state path-integral. The ground state and fluctuations around it can be obtained by this 
method. In this picture, one can obtain macroscopic quantum superposition states for attractive 
Bose gas. The coherent gap of degenerate ground states is obtained with the instanton technique. 
The phenomenon of macroscopic quantum self-trapping is also discussed. 
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INTRODUCTION 



Quantum tunneling at mesoscopic scale is one of 
the most fascinating phenomena in condensed matter 
physics. The double- well potential provides a simple 
and yet physically relevant example for studies of quan- 
tum tunneling in mesoscopic systems. Recently, there 
have been great experimental and theoretical interests in 
studying the coherent quantum tunneling between two 
Bose-Einstein condensates (BEC) in a double-well po- 
tential. It was found that the ground state changes from 
a coherent state to a Fock state as the interaction be- 
tween particles is increased in a repulsive Bose gas 0,0. 
It was numerically shown that the attractive Bose gas in 
a double-well potential has Schrodinger Cat-like ground 
states0- A similar model can be found in Refs/Q, Q 
and IJ^ 

In this paper we study quantum coherence of BEC 
in a double-well potential by mapping the two site bo- 
son model onto an anisotropic spin model in an exter- 
nal magnetic field Then the coherence properties of 
Bose system can be studied in the SU(2)-coherent-state 
path-integral representation. From the effective classical 
energy, it is easy to show that the phenomenon of macro- 
scopic quantum self-trapping (MQST) exists in both the 
repulsive and attractive interaction cases. The point sep- 
arating the coherent ground state and macroscopic quan- 
tum superposition state can be obtained analytically. 
The number fluctuation and relative phase fluctuation 
between two-well condensates are given through path- 
integral technique. Within the instanton technique, we 
obtain the tunnel splitting of degenerate ground states in 
the attractive interaction case. It is noted that the model 
presented here is general (see Eq.Q), including BEC in 
symmetric or non-symmetric double-well potential, and 
in different regimes of parameters. We emphasize that 
the quantum coherence properties of BEC depend on the 
parameters of system distinctly. 



II. PHYSICAL MODEL AND SU(2) COHERENT 
STATE 

For the Bose gas in an external potential, U (r), the 
Hamiltonian can be written in the second quantized form 
as 



dfip^ipip^Tp. (1) 



In Eq. we have used a shape-independent form for 
the atom-atom interaction with g — ^■nfi^asc/'m, where 
ttsc is the s-wave scattering length for repulsive {g > 0) 
or attractive (fl < 0) interactions. Under the two-mode 
approximation|g, -0 can be expanded as 

i> ^ (l)LbL + (l)BbR, 

where (jiL and 0^ are real, and describe the mainly left- 
well and mainly right-well populated states respectively. 
With the help of Schwinger boson representation for an- 
gular momentumi^]. 



J+ ^ Jx+ iJy = blbn, J_ = j|, = i (b[ 



Lbh- b\bR 



N 



b\bL + b\bR, J 



N fN 



(2) 



we can map the Bose model (1) to a new Hamiltonian: 



H 



(3) 
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where e = S dr<pL{R) i^-^J^"^ + U) (t)L(R), Pl(r) = 

9!dr^l(Ry i' = ~ I dr(l)LiR) ("I^V^ + <f>R^L), u = 
gjdr(j)l(j>% VLi^R) = 9 1 dr(t)l(^ii)<l>R{L)- One can easily 
estimate that |/3l(_r)| 3> |ul(_r) I- And the condition 
t' > can always be satisfied by choosing proper signs 
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(positive or negative) for real <f>L and <j)R. The Hamilto- 
nian (|3Jl describes an anisotropic spin system in an exter- 
nal magnetic field, which has been studied extensively in 
quantum tunneling in mesoscopic magnets!^. The terms 
involving N but independent of components of J are con- 
stants and can be dropped in H . The Hamiltonian Eq.Q 
is general, valid for the symmetric and non-symmetric 
wells as well. By non-symmetric well we mean the differ- 
ence in the depth and shape between left and right wells, 
which can be tuned by the external potential U (f) . 

For simplicity, we consider the symmetrical-well case, 
f^L = Pr = P' , vl = vr = V. Hence the Hamiltonian 
can be simplified a,s H ~ —2tJx + PJ^ + 2mJJ, where 
P = P' + u and t — t' — vN. In this paper, we consider 
only the case |/3| ^ |m| and t > 0. In fact, the two-mode 
approximation is invalid for large vN for the occupation 
of other modes _6j. The Hamiltonian now reads: 



(a) prj/2t-2 



(b) prj/2t-i 



H = ~2tJ^ + pjl 



(4) 



In the SU(2)-coherent-state path-integral representation, 
the Euclidean transition amplitude from an initial state 
to a final state can be written as^^ 



e-f^(^^-^')/'* = J [dn{T)] exp 



where 



Se 



dT 



N 

ih—{l — COS 6) 



rf7 



(5) 



(6) 



The first term in Eq. ((HJ is the Wess-Zumino term, and 
the effective classical energy E{9, cj)) is 



E{e,(j)) = -tN sine cos i 



PN^ 



■ cos 



(7) 



It is noted that the action © describes the 
dimensional dynamics in the Hamiltonian formulation, 
which consists of the canonical coordinates and the 
canonical momentum — ihN {l~cos9) /2. In this 
picture, {Jx) = sinScosi/), {Jy) — ^ si'O-O sm(f>, (J^) — 
§cose = {Nl - Nr)/2, where Nl and Nr are the left 
and right-well particle number respectively. Therefore, 
cos 9 = {Nl — Nii)/N corresponds to the double-well 
population imbalance, and {b]^bfi} = (J+) — Ysin6'e*"^, 
(j) corresponds to the phase difference between the double- 
well condensates. 



III. EFFECTIVE CLASSICAL ENERGY 

In this section, we will show that the effective classical 
energy and the classical equations of motion give many 
interesting results about quantum coherence properties 
of EEC in a double-well potential. It is natural to inves- 
tigate the classical orbits on the Bloch sphere, which can 




FIG. 1: Classical orbits (energy scale is tN) for different pa- 
rameters l3N/2t. cos 9 corresponds to the double-well popula- 
tion imbalance and corresponds to the phase difference be- 
tween the double-well condensates. Degenerate energy max- 
ima (or minima) imply the phenomena of MQST. Figs. 1(a) 
and (b) are for the repulsive case, Figs. 1(c) and (d) are 
for the attractive case. Two critical points correspond to 
/3N/2t = ±1. 



be described by the energy contour. In Figs. l(a)-(d), we 
plot the classical orbits for different parameters PN/2t. 
The classical orbits show the interesting phenomenon of 
self- maintained population imbalance, i.e., macroscopic 
quantum self-trapping. This phenomenon was first found 
by Smerzi et al. in both repulsive and attractive inter- 
action cases by using the canonical conjugate variables 
approximation and numerical calculation, and was ex- 
plained as a nonlinear phenomenon induced by the in- 
teraction between atoms j^. Here with the help of effec- 
tive classical energy, we study this MQST phenomenon 
in both repulsive and attractive interaction cases, and 
present the condition for MQST analytically. From the 
effective classical energy E{9, cj)) one can easily obtain the 
two dividing points PN/2t = ±1, which correspond to 
the existence of two degenerate energy maxima a,t (j) — -k 
(i.e. the phase difference between double-well BEC is tt) 
or minima at (/> = (i.e. the phase difference is 0). When 
-1 < PN/2t < 1, Fig. 1 show that MQST is forbidden 
and the particle number of each well oscillates around 
N/2. However the phenomena of MQST exist in the cases 
of PN/2t > 1 and PN/2t < -1. As a result, we conclude 
that the phenomenon of MQST is permitted when the 
interaction between atoms is strong enough to obtain de- 
generate energy maxima or minima {\P\N/2t > 1) in both 
repulsive and attractive interaction cases. Our results of 
repulsive case agree well with the results in Refs.Q and 



repi: 
g. 



Another interesting observation concerns macroscopic 
quantum superposition state in BEC with attractive in- 
teraction. One can easily show that the system has differ- 
ent energy minima for different parameters PN/2t. The 
energy minima appear at </> = and different 9 (i.e. dif- 
ferent population imbalance). The 6'-dependence of effec- 
tive classical energy E{9^ </() = 0) is plotted in Fig. 2 for 
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FIG. 2: ^-dependence of effective classical energy for different 
parameters PN /2t. 



different parameter f3N/2t. When (3N/2t > -1, there 
is only one energy minimum at 9 = 7r/2; while when 
(3N/2t < — 1, there are two degenerate energy minima at 
61 = 00 and 6* = TT - 6^0, where sin 6*0 = 2t/\l3\N. The for- 
mer case favors a coherent ground state or a Fock ground 
state0, the latter case favors a Schrodinger cat-like 
ground state which is the superposition of two SU(2) co- 
herent states 16*0,0) and [tt — 6*0,0). This kind of defini- 
tion describes the same Schrodinger cat-like state as that 
in Ref. 3] , where the result was obtained by numerical 
calculation. Here the dividing point separating the co- 
herent ground state and the Schrodinger cat-like state 
is analytically obtained in a simple and clear approach. 
Moreover, the tunnel splitting can be obtained by apply- 
ing the instanton technique in the SU(2)-coherent-state 
path-integral representation, as shown in the next sec- 
tion. 



IV. COHERENT QUANTUM TUNNELING 

As shown in Fig. 2, there will be two degenerate en- 
ergy minima when (3N/2t < — 1 for the attractive BEC. 
Now the system in question performs coherent quantum 
tunneling (i.e., quantum coherence or coherent super- 
position) between two degenerate energy minima. The 
tunneling removes the degeneracy of the original ground 
states, and the true ground state (i.e., Schrodinger cat- 
like state) is a superposition of the previous ground 
states. Tunneling between neighboring degenerate vacua 
can be described by the instanton configuration and leads 
to a level splitting of the ground states. Here we evaluate 
the tunnel splitting of two degenerate ground states by 
applying the instanton technique^]. 

After adding some constants, we can rewrite the effec- 
tive energy as 



(sin 6* — sin 6*0)^ + tN sin (1 — cos( 



where sin 6*0 = 2t/\P\N. From SSE{0,(f>) = 0, we obtain 
the instanton solution as 



cos 9 
sin (j) 



cos 9o tanh (w^r) , 
i cot^ ^0 sech^ (wbr) 
[1 + cot^ 6*0 sech^ (wfcr)] 



1/2 ' 



(9) 



corresponding the transition from 9 = 9o to 9 = tt — 
9q, where ojb = N\P\cos9o/2h. The associated classical 
action is found to be 



Sc 



N 



1 f l + cos9o 
cos 9n -\ — in 

2 V 1 - cos 6*0 



(10) 



and the final result of tunnel splitting is 



HA = 8 



271 1 sin 6*0 



1 — cos 6*0 
1 + cos 6*0 



^ cos &0 



-Sal 



(11) 

It is noted that the tunnel splitting is obtained with 
the help of instanton technique in the SU(2)-coherent- 
state path-integral representation, which is semiclassi- 
cal in nature, i.e., valid for large N. Therefore, one 
should analyze the validity of the semiclassical approx- 
imation. The semiclassical approximation is valid only 
when the energy splitting is far less than the energy bar- 
rier A^^|/3| (1 — sm9o)^ /4 and the energy of zero point 
oscillation -/V^|/3| cos^ 9o/A, which indicates that the clas- 
sical action Sd ^ 1. From Eq. (|10|l one can easily see 
Sci 3> 1 when 9q — 7r/2 — 0.3 for typical particle num- 
ber of attractive BEC N = 1000. Then the semiclas- 
sical approximation should be already rather good for 
0<9o< 9i 

Because of coherent quantum tunneling, the two 
degenerate energy minima ground state |0o,O) and 
\n — 9o,0), which correspond to SU(2) coherent states 
with population imbalance ± cos 9o and phase differ- 
ence between condensates in the two wells, split to 
the two parity-different Schrodinger cat states |±) = 
(|6lo, 0) ± Itt - 6*0, 0)) /V2 with energy splitting hA. One 
can show that the energy splitting is extremely small 
when ^0 is away from tt/2. 



V. NUMBER FLUCTUATIONS AND RELATIVE 
PHASE FLUCTUATIONS 

In this section we consider the fluctuations around 
the ground state. We shall evaluate the fluctuations for 
the parameters 9 and (j), which correspond to the rel- 
ative number fluctuations and the phase fluctuations. 
Rewriting the parameters as 9{t) — 9{t) -\- 9i{t) and 
(j>{T) — <j){T) + (/)i(t), one obtains the Euclidean action as 



Se{9,(I^)^Sci + ^S^S, 



(12) 



(8) 



where Sd is the classical action which satisfy SSd — 
and 
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-6^S = 
2 



in — (sin( 

2 dr ^ 



1 

2 _ 

1 rt 



ih — COS 9 (b O^dr 
2 ^ 



(13) 



{EeeOl + 2Eg^ei4)i + i^^^^?) dr. 



In the above equation, Egg ~ d^E/dO^, Eg^ = 
d^E/dOdcf) and E^^ = d'^E/d4>^, which are evaluated 
at the classical path. Under the condition that E^^ > 0, 
the Gaussian integration can be performed over (pi , then 
the effective action for 9i is found to be 



= j i^Adl + Be,e, + cefj dr. (i4) 

Using Eq. ((HJ, we have 

— tN sin 6 cos ^ > 0, 
h^N sine 

A 



8t cos (j) 
B ^ 0, 
^ tN cos (f) 

The effective action reads: 



sm ( 



fi^Nsme ■ 



8t cos (j) 
tN cos PN^ 
2 sin 9 ^ 



dr. (15) 



In the case of PN/2t > — 1, the classical ground state 
is ^ = 7r/2, = 0. Now we study the 9 fluctuation near 

the classical ground state, A = and C = ^ + = 
f (2i + /37V). The effective action is 



8t 



N 



{2t + i3N)el 



dr. 



(16) 



The motion of the fluctuation iV^i is approximately 
a harmonic oscillator with mass 



■^-^ and frequency 
?:!^21+J^I11_ Jts characteristic length crjve is determined 



by (Jn0 = 
one well is 



N 



2t 



2t+PN ■ 



Hence the number fluctuation in 



In the case of j3N/2t < — 1, the classical ground state 
is ^ = 6*0 or TT - 6*0, ^ = 0,where sin 6*0 = 2t/\(3\N. How- 
ever, due to quantum tunneling between the two classi- 
cal ground states, the true ground state is an even-parity 
Schrodinger cat state. Inspecting the fluctuation near the 
classical ground state, we find that A = ^"^g"/"^" = ^ 



and C = .ri^ 

2 sin 9 

five action is 



I3N^ 
4 



4|/3| 



cos^ 6*0. The effec- 



^COS^.0^? 



dr. (18) 



The motion of the fluctuation A^^i is approximately 
a harmonic oscillator with mass .^i " ^„ and frequency 



,3|jVcos 00 
h 

by o-irg = 
well is 



2|/3|Af2 

Its characteristic length a^e is determined 



2^ 



Hence the number fluctuation in one 



A {Nl - Nr) = Nsin9oA9 = sin 6*0 



2N 
cos 6*0 ' 



(19) 



which is different from the case of (3N/2t > —1. This 
result is not the true number fluctuation around the 
Schrodinger cat state. It indicates the fluctuation around 
the classical ground state. The true number fluctuation 
is of the order of N'^, referred to Ref. nj] as a superfrag- 
mented state. 

Finally we discuss the relative phase fluctuation. After 
performing the Gaussian integration over 9i , the problem 
reduces to the one-dimensional path integral with the 
effective action: 



with 



A' 

C ■ 



A'(j>l + C'^l)dT, 



(20) 



S^TV^sin^ 



{Egg- cot 9Eg)' 
El 





In the case of f3N/2t > — 1, the classical ground state 
is ^ = 7r/2, ^ = 0. Now we inspect the (j) fluctuation near 
the classical ground state, A' ~ ^ ^ 
The effective action is 



4{2t+f3N) 



and C = if. 



A(iVi-7V^)= (7Vsin|) 



A9 




(17) 



One can see that there is a singular point in f]N/2t — 
— 1, which indicates a dividing behavior between coherent 
state-like and Schrodinger cat state-like ground state. 



h^N 



4 {2t + PN) 



^1 + 



The motion of the fluctuation 
harmonic oscillator with mass 



tN 



dr. 



(21) 



bi is approximately a 
^ and frequency 



2(2t+/3Af) 
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^liZl^fim^ Its characteristic length is determined 



by o'0 = • Hence the relative phase fluctua- 



tion is 




(22) 



In the case of f3N/2t < — 1, the classical ground state 
is ^ = 6*0 or TT - So, ^ = 0, where sin6lo 2t/\(3\N. 
Inspecting the fluctuation near the classical ground state, 
we find that A' = 4^f^",'/g^ and C = \l3\N^ sin^ Oq. The 
effective action is 



sin^ 00 
A\(3\ cos2 6*0 



+ \/3\N^ sin^ 9o<j>l 



dr. 



(23) 

The motion of the fluctuation 4>i is approximately a 
harmonic oscillator with mass 2\i3 \ cos^eo frequency 
2|/3|jVcos 00 ^ characteristic length is determined by 
2 _ — coseo — Hence the relative phase fluctuation is 



a 



2N sin^ So ' 



: 



/ COS Ua 
2iVsin2 9o' 



(24) 



VI. CONCLUSION 

In summary, we study the quantum coherence phenom- 
ena in a double- well BEC based on SU(2)-coherent-state 
path-integral. By this method we analytically study the 
MQST phenomenon, the ground states of this system, 
and the existence of macroscopic quantum superposition 
states. We find that MQST will happen in both repulsive 
and attractive interaction cases, and analytically obtain 
the dividing points l3N/2t = ±1. When l3N/2t > -1, 
both repulsive and attractive Bose gases favor a coherent 



or a squeezed ground state, which is the coherent state for 
non-interaction case, the relative-number-squzeed state 
for rcplusive case, and the relative-phase-squzeed state 
for the attractive case (see Eqs.(17) and (22)) respec- 
tively. However when /3N/2t < — 1, attractive Bose gases 
favor a macroscopic quantum superposition state. The 
relative number fluctuation and relative phase fluctua- 
tion between two-well condensates are obtained through 
path-integral technique. For the attractive interaction 
case, the coherent gap of degenerate ground states is ob- 
tained analytically with the help of the instanton tech- 
nique. 

It is noted that all the discussions in this paper does 
not impliy this kind of macroscopic quantum superposi- 
tion states can be easily created. To inspect the decoher- 
ence, one has to induce the interaction between thermal 
could and condensate. Dalvit et al. have studied the de- 
coherence of a similar model in Ref > which is beyond 
our two-mode approximation and semiclassical approach. 
However, as an rough estimate, the high barrier and small 
splitting gap imply high rate of decoherence. In fact, to 
the best of our knowledge, there is no experimental ev- 
idence of BEC cats up to now, although there are some 
proposals to create them|^[l^[T?]. Ref. 12] pointed out 
that such macroscopic quantum superposition states are 
extremely fragile to decoherence, and suggested that the 
strategy of trap engineering and symmetrization of the 
environment should be able to deal with that issue. The 
theoretical calculations performed in this paper can be 
extended to the Bose gas in the non-symmetric double- 
well potential, and an effectively two-component spinor 
condensate. Work along this line is still in progress. 
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